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ABSTRACT 

The Sun will eventually lose about half of its current mass nonlinearly over 
several phases of post-main sequence evolution. This mass loss will cause any 
surviving orbiting body to increase its semimajor axis and perhaps vary its ec- 
centricity. Here, we use a range of Solar models spanning plausible evolutionary 
sequences and assume isotropic mass loss to assess the possibility of escape from 
the Solar System. We find that the critical semimajor axis in the Solar System 
within which an orbiting body is guaranteed to remain bound to the dying Sun 
due to perturbations from stellar mass loss alone is « 10 3 AU - 10 4 AU. The fate 
of objects near or beyond this critical semimajor axis, such as the Oort Cloud, 
outer scattered disc and specific bodies such as Sedna, will significantly depend 
on their locations along their orbits when the Sun turns off of the main sequence. 
These results are applicable to any exoplanetary system containing a single star 
with a mass, metallicity and age which are approximately equal to the Sun's, 
and suggest that few extrasolar Oort Clouds could survive post-main sequence 
evolution intact. 

Key words: planet-star interactions, planets and satellites: dynamical evolution 
and stability, stars: evolution, stars: AGB and post-AGB, Oort Cloud, minor 
planets, asteroids:general 



1 INTRODUCTION 

The fate of our Solar System is of intrinsic human inter- 
est. From the present day, the Sun will continue to exist in 
a main sequence phase for at least 6 Gyr longer, and then 
experience post-main sequence phases before ending its 
life as a white dwarf. The post-main sequence phases will 
be violent: the Sun's radius and luminosity will likely vary 
by several orders of magnitude while the Sun ejects ap- 
proximately half of its present-day mass through intense 
winds. The effect on orbiting bodies, including planets, 
asteroids, comets and dust, may be catastrophic. 

The future evolution of the Solar System between 
now and the end of the Sun's main sequence will be 
largely unaffected by the Solar mass loss. The Sun's cur- 
rent mass loss rat e lies in the range ~ 10~ 14 M(T)/yr - 
lO- 13 M /yr (e.g. ICohenl bulll: IPitieva fc Pitievl [201ll V 



As demonstrated by IVeras et al.l (|2011l ). all orbiting ob- 
jects would evolve adiabatically due to this small mass 
loss. Thus, assuming the mass loss rate remains con- 
stant by the end of the main sequence, their eccentric- 
ities will remain unchanged, and their semimajor axes 
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would increase by at most ~ 0.055%. The Sun's main se- 
quence mas s loss rate, has, however, varied o ver time. 
IWood et all l|2002h and IZendeias et ail (|2010t ) suggest 
that this mass loss rate was orders of magnitude greater 
when the Sun began life on the main sequence, and is 
a monotonically decreasing function of time. Therefore, 
0.055% represents an upper bound. Regardless, this po- 
tentially slight orbital expansion, uniformly applied to all 
orbiting objects, is not predicted to change the dynamics 
of the Solar System. 

Instead, the primary driver of dynamical change in 
the Solar System during the Sun's main sequence will 
arise from the mutual secular perturbations of the orbit- 
ing bodies. Ikholshevnikov fc Kuznetsovl l|2007l ) provide a 
comprehensive review of the investigations up until the 
year 2007 which have contributed to our understanding 
of this evolution. These studies describe orbital evolution 
beyond 10 4 yrs from now in a primarily qualitative man- 
ner because numerical integrations typically cannot re- 
tain the orbital phase information of the terrestrial plan- 
ets over the Sun's main sequence lifetime. Furthe r, Mer- 
cury might suffer a close encounter with Venus (jLaskar 
bOOSt ). and the orbital motion of the other inner planets 
is predicted to be chaotic. The outer planets, however, 



© 2012 RAS 



2 Veras & Wyatt 



are predicted to survive and harbor quasi-periodic mo- 
tion. Since this review, additional studies have primarily 
foc used on the complex evolution of the terrestrial plan- 
ets; [Laskarjjtjil] (|201 lh helps affirm the intractability of 
terrestrial orbital phases o ver long times. O ther studies 
l|Batygin fc Laughlinll2008l and lLaskadl^OOSh corroborate 
previous results about how the outer Solar System plan- 
ets (Jupiter, Saturn, Uranus, and Neptune) should re- 
main stable in near-circular (e < O.f) orbits, even though 
the orbits themse lves may be chaotic !|Havesl l2007, 2008; 
lHaves et al1l2010h . 

This uncertainty of a planet's orbital architecture at 
the beginning of the Sun's post main sequence phases 
renders difficult the determination of the fate of individ- 
ual presently-known bodies. Regardless, several studies 
have explored the effect of the Sun's post main-sequence 
evolution on the planets. The outcome for the Earth is 
controversial because investigators differ on whether the 
Earth's expanding orbit will "outrun" the Sun's expand- 
ing envelope, a nd on how to model the resulti ng tidal 
interactions fsee lSchroder fc Connon Sm ith 2008 and ref- 
erences therein - ) . iDuncan fc Lissauerl l|l998h consider the 
stability and orbits of Jupiter, Saturn, Uranus and Nep- 
tune with mass scalin gs and a post main-seq uence mass 
loss prescription from lSackmann et al.l (|l993). They find 
that this prescription yields "no large growth in plane- 
tary eccentricities", but mention the possibility of Oort 
Cloud comets being ejected during periods of brief, rapid 
mass loss. Complimentary studies of other planetary sys- 
tems have explored in more detail ho w Oort Cloud comet- 
analogues may become unbound l|Alcock et alj 1 19861 ; 
iParriott fc Alcocklll998r i. 

Our study is devoted to exploring the prospects for 
dynamical ejection from the Solar System by using real- 
istic multi-pha sic nonlinear mass l oss prescriptions from 
the SSE code (|Hurlev et all 1 2000) and by treating iso- 
lated objects beyond the influence of any potentially sur- 
viving planets . We do so by using the analytical mass loss 
techniques in IVeras et al.l (|201lh , which yield a critical 
semimajor axis, a cr i t , within which an orbiting object is 
guaranteed to remain bound. In Section 2 we present the 
solar evolutionary models used, highlighting those which 
are later applied to N-body simulations that help reaffirm 
the validity of a cr it- In Section 3, we define a cr it, plot it 
as a function of Solar evolution model, and explore the 
properties of some representative orbits. In Section 4, we 
explore the implications for the scattered disc and the 
Oort Cloud. We discuss the results in Section 5 and con- 
clude in Section 6. 



2 SOLAR EVOLUTION MODELS 



We utilize the SSE code (|Hurlev et al.ll2000h . which yields 
a complete stellar evolutionary track of a star for a given 
initial stellar mass, metallicity and mass loss prescrip- 
tion. We denote the Sun's mass as M*(t), and define a 
Solar Mass = 1M© to be the current value of the Sun's 
mass. We assume that the Sun's current metallicity is 
[Fe/H]© = 0.02. We also assume that the Sun has al- 
ready been evolving on the main sequence for exactly 4.6 
Gyr. 
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Figure 1. Potential Solar evolutionary tracks. The upper 
panel displays at least a part of all phases from the RGB 
phase to the thermally pulsing AGB phase for 7 representa- 
tive values of tj© (the nearly-flat RGB curves before 12.290 
Gyr are not shown). The lower panel zooms-in on and isolates 
the thermally pulsing AGB phase for 77© = 0.2 (solid line), 
r?© = 0.3 (dotted line) and t)q = 0.4 (dashed line). Time is 
measured from the beginning of the Sun's main sequence. The 
plot demonstrates how the choice of r/© can drastically alter 
the percentage of Solar mass lost and Solar mass loss rate per 
phase. 



The remaining, and unknown, constraint to include 
is the Sun's post-main sequence mass loss. On the 
Asymptotic Gian t Branch (AGB), we use t he mass-loss 
prescription from IVassiliadis fc Wood! (| 19931 ): 



log 



dAL 
dt 



= -11.4+0.0125 [P - lOOmax (A/* - 2.5, 0.0)] (1) 



such that dh'h/dt is computed in M©/yr and where 

logP = min(3.3, -2.07 - 0.9 log M* + 1.94 log R+) , (2) 

where P is computed in years and R+ is the radius of the 
Sun. 

On the Red Giant Branch (RGB), mass loss is often 
mode led to have the following "R eimers Law" functional 
form jKudritzki fc Reimerdll978r ): 



dM * - r, (4 x HT 13 ) L ^)Mt)M e 



(3) 



dt y ' M t yr 

where represents the Solar luminosity and r\ is an im 
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portant dimensionless coefficient whose magnitude helps 
determine the sequenc e of stellar evolutionary phases. 
ISackmann et alj (|l993h compute evolutionary models for 
the Sun, and state that 0.4 < 77 < 0.8 represent "reason- 
able mass-loss rates". This c onforms with the typically- 
used value of 77 = 0.5 (e.g. Kudritzki fc Reimerj Il978l ; 
iHurlev et all l200d ; iBonsor fc Wvattll2010l). An updated 
version of the Reimers Law (|Schroder fc Cunt j 120051 ') 
contains two new multiplicative factors: 



3 THE CRITICAL SEMIMAJOR AXIS 
3.1 Equations 

A star which is expunging mass beyond an orbiting object 
will cause the latter's orbit to expand. If the mass loss is 
large enough during one orbital period, then deformation 
of the orbit can cause the object to escape the system. We 
quantify th ese claims by defin ing a dimensionless "mass 
loss index" (Verasetal 
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where <7* and T e g are the star's surface gravity and effec- 
tive temperature, and <?© is the current value of the Solar 
surface gravity. The product of both t hese additional fac- 
tors is not far from unity. According to lSchroder fc Cunt j 
(2005), who used the RGBs and Horizontal Branches of 2 
globular clusters with very different metallicities for cali- 
bration, their effective value of 77 in the original Reimers 
law would be about 0.5. 

All these considerations dissuade us from selecting 
a single value of 77© and instead lead us to consider the 
entire range 0.2 ^ t]q ^ 0.8; we sample 770 values in this 
range in increments of 0.01. Representative evolutionary 
sequences in this range are displayed in the upper panel 
of Fig. [T] In all cases, the Sun undergoes a > 0.5 Gyr- 
long transition towards the foot of the RGB after leaving 
the main sequence. The Sun moves on to the RGB before 
enduring a period of core helium burning a nd subsequent 
evolut ion on the early AGB, as defined by IHurlev et al.l 
(120001 ). Then, for r) Q ^ 0.74, the Sun will enter a ther- 
mally pulsing AGB phase before becoming a white dwarf. 
This intermediate phase is bypassed for 0.75 ^ r/Q ^ 0.80. 

The different evolutionary sequences demonstrate 
that the greatest amount of mass lost and the greatest 
mass loss rate may occur during the RGB, Early AGB, 
or Thermally Pulsing AGB phases. In all cases the to- 
tal percent of the Sun's current mass which is lost is 
~ 46.5%-49.0%. The value of 77 determines how much en- 
velope mass is lost on the RGB, and hence, the amount 
of envelope mass which remains at the start of the AGB. 
For 770 < 0.5, most of this disposable mass is reserved 
for the AGB, whereas for 770 > 0.5, most of this mass 
is lost during the RGB phase. When 77 is large enough 
(770 ^ 0.75), then enough of the envelope mass is lost dur- 
ing the RGB to cause the Sun to bypass the Thermally 
Pulsing AGB phase. In the extreme case of 770 ^ 0.94 
(not modeled here), the Sun would bypass Helium burn- 
ing and the AGB phases entirely, and would transition to 
a pure Helium white dwarf directly from the RGB phase. 

The lower panel zooms in and isolates the thermally 
pulsing AGB phase for 770 = 0.2 — 0.4, where the mass 
loss rate is greatest. Notice that the timescale for this 
mass loss is ~ 0.1 — 0.3 Myr, orders of magnitude shorter 
than the mass loss on the RGB. 



mass loss timescale a 
v|> = = 

orbital timescale n/i 

= 1 ( - \( a ^ ( M V 1 (S) 
2tt \lM Q /yr J V1AU7 \IM Q ) ' y ' 

where a represents the Solar mass loss rate, a and n 
represent the orbiting object's semimajor axis and mean 
motion, and fi = M* + M s , where M 3 represents the 
mass of the object. We can assume /j « M* because 
every known object which orbits the Sun has a mass 
that is < 0.1% of IMq. However, observations cannot 
yet exclude the possibility of the existence of a Jovian 
planet orb iting at a distance of at least several thou- 
sand AU l|lorid [201 ll ). Therefore, "undiscovered" plan- 
ets beyond ~ 10 4 AU might be massive enough to make 
a nonnegligible contribution to fi. Even if such a planet 
formed within the orbit of Neptune, gravitational scatter- 
ing could have pla ced the planet in its present location 
(jVeras et al J 120090 . Alternatively, the planet could have 
previously been free-floating and subsequently captured 
by the Solar System. 

Dynamical evolution critically depends on fy. When 
$ < 1, the orbiting object evolves "adiabatically" such 
that a increases but its eccentricity, e, remains constant. 
Examples of objects with $ < 1 include all 8 of the 
inner and outer planets, both currently (see Section 1) 
and during any plausible post-main sequence Solar evo- 
lution phase. Alternatively, when ^ 3> 1, in the "run- 
away" regime, a continues to increa se but now e may 
achieve any value from zero to unity. IVeras et all (I2OIII ) 
demonstrate that the bifurcation point beyond which the 
motion is no longer adiabatic is not sharp and occurs 
at ^bif ~ 0.1 — 1.0. In fact, if is instead defined as 
aT/n, where T is the orbital period of the orbiting ob- 
ject, then this definition differs from Eq. (JSJ by a factor 
of l/(27r) « 0.16. Therefore, both definitions effectively 
bound the approximate location of the bifurcation point 
when ^ = 1, and we consider both possibilities in our 
analysis. 

The Sun loses mass nonlinearly, as showcased in 
Fig. [[] One can ap ply the analytical linear results from 
IVeras et alJ l|201lf ) to a nonlinear mass loss profile by 
partitioning the profile into N smaller approximately lin- 
ear stages, with mass loss rates of ai, 0:2, ...a;...Qjv. If all 
these segments occur in the adiabatic regime, then we can 
relate the semimajor axis and eccentricity at the start of 
the first segment to these same elements at the end of 
segment i with: 

a» = a ( — ] sao/3j~\ (6) 
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eo, 



(7) 



3.3 N-body Simulations 



Because the orbiting object's eccentricity does not change 
in this regime, the body will remain bound in the absence 
of additional perturbations. 

All orbiting bodies will be expanding their semima- 
jor axis at the same rate. The mass loss index at the end 
of segment i will be: 



"5; — net. 



f aoMo\ ' 



(8) 



where k = 1 or 2ir, depending on whether <& is defined 
with respect to the mean motion (Eq. [5]) or the orbital 
period, k could also be treated generally as a tunable pa- 
rameter which may better pinpoint the location of the 
bifurcation point in specific cases. The first instance at 
which the orbiting object may escape is when the evo- 
lution becomes non-adiabatic and hence the eccentricity 
can vary. This occurs when ~ 1. Therefore, the min- 
imum percent of no which must be retained in order to 
guarantee that the orbiting object will remain bound in 
a given segment i is: 



Pcrit; 
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(9) 



Hence, the critical oo, termed a c rit, at which an orbit- 
ing object will remain in the adiabatic regime, and thus 
bound to a 1 Mq star, is 



1AU 



k 3 mm 



lM Q /yr 



(10) 



where the minimum is taken over all segments i. 



3.2 Mass Loss Profile Partition 

The accuracy of a cr i t depends on how finely the Solar 
mass loss profile is partitioned into stages. Given the fo- 
cused nature of our study, we simply adopt the most accu- 
rate option by treating each SSE timestep as a segment. 
These segments also include transitions between evolu- 
tionary phases, for which a mass loss and timescale is 
computed by SSE in every case. We choose output param- 
eters which yield 3333 outputs total per simulation for 
the RGB, core helium burning, early AGB and thermally 
pulsing AGB phases, and at least 20 outputs for each 
of the other phases of evolution. Our choices are moti- 
vated by the observation that the segment which satisfies 
Eq. (|10|l always arises from the RGB or an AGB phase. 
The minimum SSE timesteps are « 10 3 yr long, which is 
shorter than the orbital period of the objects we consider 
here. Equation (|10|) is valid for any relation between SSE 
timesteps and orbital period because, in the adiabatic 
regime, the semimajor axis and eccentricity evolution is 
independent of the location of the object along its orbit. 
We apply this output parametrization uniformly across 
the t]q values sampled. 



We supplement the computation of a cr it (from the stel- 
lar evolution profiles alone) with N-body simulations 
of orbiting objects. These simulations help affirm that 
objects with a < a cr i t do remain bound, and pro- 
vide a rough characterization of the motion of ob- 
jects which are susceptible to escape. We select the 
7 Solar evolutionary pathways corresponding to r/Q = 
[0.2,0.3,0.4,0.5,0.6,0.7,0.8] for these simulations, and 
use a second-order mixed-variable symplectic integrator 
from a modified vers ion of the Mercury integration pack- 
age |Chambers| [l999'). All orbiting objects are treated as 
test particles for computational ease and so that they 
could be integrated simultaneously. Because n is dom- 
inated by M^, an isolated giant planet would react to 
stellar mass loss nearly equivalently to an isolated test 
particle on the same orbit. After each Mercury timestep, 
we linearly interpolate the value of the Sun's mass from 
the SSE output. This approximation is sufficient for nu- 
merical integrations which don't feature close encounters. 

Simulating orbital evolution throughout the remain- 
ing lifetime of the Solar System with numerical inte- 
grations and without utilizing secular approximations is 
computationally unfeasible and, for this project, unnec- 
essary. Our focus is on the behavior of objects when they 
are most susceptible to escape. Therefore, having already 
identified when the greatest mass loss occurs, we simu- 
late the entire thermally pulsating AGB phase for the 
r/Q — [0.2,0.3,0.4] simulations (see Fig. [TJ. The duration 
of these simulations is between 2.7 x 10 s yr and 5.7 x 10 5 
yr. For r/Q = [0.5,0.6,0.7,0.8], we simulate Solar mass 
evolution from when M* = 0. 98000 Mq, which occurs 
during the RGB, until the beginning of the white dwarf 
phase (see Fig. [TJ. For these simulations, modeling the 
phases after the RGB may be important because despite 
their overall low mass loss rates, the AGB phases could 
feature short bursts of high mass loss. Also, /?,: is lower 
for the phases beyond the RGB. The duration of these 
simulations is between 1.64 x 10 8 yr - 1.88 x 10 s yr. 

We choose the initial conditions for the orbiting ob- 
jects as follows. The orbital evolution of an orbiting ob- 
ject which is susceptible to esca pe is highly depen dent 
on e and its true anomaly, / fsee lVeras et al.ll201l| y For 
linear mass loss and 9 3> 1, the object's eccentricity 
is likely to initially increase. However, the eccentricity 
must initially decrease if / cr it ^ / ^ (360° — /crit), where 
/ crit = 180° - (1/2) cos" 1 (7/9) « 160°, and we wish to 
sufficiently sample this behavior as well. These consider- 
ations motivate us to choose /o in increments of 30° from 
0° - 330°, supplemented by f = [160°, 170°, 190°, 200°]. 
In order to sample a wide range of eccentricties, we choose 
e = [0.0,0.2,0.4,0.6,0.8,0.95]. Because application of 
Eq. (HUJ) demonstrates that 10 3 AU < a cr it < 10 4 AU, we 
choose a /AU= [10 2 ,10 3 ,3 x 10 3 ,5 x 10 3 , 10 4 , 10 s ]. The 
smallest of these values yields a period of 414 years for 
M* = 0.5Mq. Therefore, we choose a timestep less than 
l/50th of this value, 1 x 10 4 days, for all our numeri- 
cal simulations. Because this timestep might not sample 
the pericenter sufficiently for highly eccentric orbits, we 
perform an identical set of simulations with a timestep 
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The Solar System's Stability Boundary 
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Figure 2. The Solar System's critical semimajor axis, as a function of the Solar Reimers Mass Loss coefficient, tjq. The curve is 
derived from Eq. (HOP . In the shaded pink region above the curve, orbiting bodies may escape the dying Sun on hyperbolic orbits. 
In the white region below the curve, bodies are guaranteed to remain bound on elliptical orbits. The colour of the line segments 
indicates the Solar phase in which orbiting bodies are most susceptible to escape: red = RGB, green = early AGB and blue = 
thermally pulsing AGB. Checkmarks and crosses provide information about the results of N-body integrations (see Section 13.311 . 
Checkmarks indicate that every orbiting object in every N-body simulation remained bound. Each cross represents a different 
initial eccentricity value (out of the set {0.0, 0.2, 0.4, 0.6, 0.8, 0.95}) of at least one body which escaped. An alternate curve, with 
k = 2ir, would have the same form but with values which are (27r) 2 / 3 PS 3.4 times lower than the ones shown here, yielding a more 
conservative estimate. This plot demonstrates that Solar System bodies must satisfy a < 10 3 AU when the Sun turns off of the 
main sequence in order to be guaranteed protection from escape. 



which is one order of magnitude lower as a check on the 
results. 



3.4 Visualizing the Critical Semimajor Axis 

We plot a cr i t from Eq. (|10[) assuming k — 1 as a function 
of 770 in Fig. [2] The curve illustrates that for 770 ^ 0.5, 
a crit ~ 10 3 AU, and for t] 3 = 0.7, a crit ~ 10 4 AU. Other 
values of a cr i t lie in-between these two extremes. The 
curve is coloured according to the stellar evolutionary 
phase at which an orbiting object is most susceptible to 
ejection, which depends on a combination of the mass 
loss rate and the mass remaining in the Sun. Three dif- 
ferent phases (red = RGB, green = early AGB, blue = 
thermally pulsing AGB) are represented on the curve, 
demonstrating the sensitivity of a cr it to the stellar model 
chosen. When the RGB evolution dominates the orbit- 
ing object's motion, the critical semimajor axis is higher 
than in the AGB-dominated cases because RGB evolu- 
tion takes place earlier in the Sun's post-main sequence 
life, when the Sun has more mass. Any objects in the pink 



shaded region above the curve might escape; those in the 
white region below the curve are guaranteed to remain 
bound. 

The overlaid symbols indicate the results of the N- 
body simulations. Checkmarks indicate that every or- 
biting object from every {eo,/o} bin sampled remained 
bound (e < 1) for the duration of the simulation. Crosses 
indicate that at least one orbiting object become un- 
bound, either through the orbit becoming hyperbolic 
(e > 1) or through the body achieving a > 10 6 AU. 
The number of crosses indicate the number of eo bins in 
which at least one orbiting object became unbound. The 
N-body simulations help affirm that all of the objects 
below the curve remained bound. This also holds true if 
k — 2n. In that case, the curve would keep the same form 
but be lowered by a factor of (2ty) 2 ^ 3 ~ 3.4. This value 
of k would provide a more conservative estimate for the 
Solar System's critical semimajor axis, and would still 
ensure that any known surviving planets and belts would 
remain bound. Objects with high values of eo are usually, 
but not always, prone to escape; additional cross symbols 
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suggest that objects escaped at progressively lower values 
of eo- The number of crosses displayed do not show uni- 
form patterns primarily because of the complex depen- 
dencies of the evolution in the adiabatic/runaway tran- 
sitional regime on a, e and / (see Section 4.2). Also, the 
parameter space sampled is limited and the N-body sim- 
ulations do not model the entire pre-white dwarf lifetime 
of the Sun. 

Not shown in Fig. [5] are the results from the ao — 10 2 
AU and ao = 10 AU simulations. In the former, every 
orbiting object remains bound (corresponding to a row 
of checkmarks). In the latter, at least one orbiting object 
escapes for each value of t/q integrated (corresponding to 
a row of different numbers of crosses). 



3.5 Properties of the Critical Stage 

Here we take a closer look at the critical stage which 
yields a cr it through Figs. 1 3161 Figures 1 3141 demonstrate 
that the critical stage occurs nearly at the end of the 
Sun's post-main sequence life, when the Sun will have 
already lost the majority of its disposable mass. The solid 
black curve in Fig. [3] is a poor proxy of the a cr it curve 
from Fig. [2] Both Figs. [3] and 0] suggest that 770 = 0.7 
represents a bifurcation point in Solar evolution, which 
is approximately where the thermally pulsing AGB stage 
becomes transitory. The critical time from the start of 
the Sun's main sequence will be ~ 12.465 ±0.005 Gyr for 
77© < 0.7, and 12.326±0.0002 Gyr for 770 > 0.7. These are 
remarkably well-constrained values relative to the curves 
in Fig. H 

Now consider a; at the critical stage. Figure [5] shows 
that this mass loss rate varies between pa 2 x lO _7 M0/yr 
and pa 5 x lO _6 M0/yr. Also, this rate closely mirrors the 
maximum rate achieved throughout the simulation, given 
by the dashed, brown curve. Therefore, by itself, max(a?;) 
is an excellent indicator of the stability prospects of an or- 
biting body. Also, this trend suggests that short violent 
outbursts from planet-hosting stars can jeopardize the 
survival of orbiting bodies. The Solar System is subject 
to the same danger even though the Sun is not expected 
to experience a violent outburst with a > 5 x 1O -6 M0 /yr 
from these models. If we are underestimating the variabil- 
ity of the Sun after it turns off of the main sequence, then 
orbiting bodies can be in greater danger of escape than 
we currently expect. 

If such variability is on a timescale which is shorter 
than the timesteps used in these simulations, then a cr it 
is being overestimated. Figure |S] illustrates the timestep 
from which the mass loss rate at the critical stage was 
computed. In all cases, this timestep is between approxi- 
mately 980 yr and 5860 yr. The Sun's current rotational 
period is less than one month; during a single rotation, 
prominences may erupt and sunspots may appear, among 
other phenomena. The Sun's (unknown) post-main se- 
quence variability will likely be more violent, but on a 
similar timescale. If so, then short bursts of intense mass 
loss could be important but are not modeled by the com- 
paratively long SSE timesteps. 
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Figure 3. The percent of \Mq remaining in the Sun at the 
moment when orbiting bodies are most likely to escape (solid 
black curve) and at the moment the Sun becomes a white 
dwarf (dashed brown curve). Orbiting bodies are most likely 
to escape the Solar System after the Sun has lost at least 
O.35M . 
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Figure 4. The time from the start of the Sun's main sequence 
that orbiting bodies are most likely to escape the Solar System. 
Despite the complex variance in Solar evolution as a function 
of 770, this critical time has a relatively smooth, piecewise 
dependence on Mq. 

3.6 Orbital Properties 

3. 6. 1 Overview 

Here we analyze some representative orbits from the N- 
body simulations. Because the timescale for mass loss 
from the bottom panel of Fig. [1] is orders of magnitude 
shorter than the top panel, the difference in orbital prop- 
erties is pronounced. Specifically, for ao > 10 3 AU and 
rjQ < 0.5, all eccentricity variations will often occur dur- 
ing a single orbit. Alternatively, for rjQ ^ 0.5, these vari- 
ations occur over many orbits. We consider one case from 
the former category and two cases from the latter in de- 
tail. Additionally, we consider the eccentricity evolution 
for one case in the runaway regime, with ao = 10 5 AU and 
rjQ = 0.2, and demonstrate here that objects near their 
apocenters do remain bound in this CELS6, clS predicted by 
the linear theory. 
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Mass Loss Rate Between SSETimesteps 
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Figure 5. The mass loss rate at the moment when orbiting 
bodies are most likely to escape (solid black curve) and when 
this value takes on a maximum (dashed brown curve). The 
excellent agreement between the two curves suggests that mass 
loss rate alone, regardless of the amount of mass remaining in 
the Sun, is the primary indicator of when an orbiting body is 
no longer safe from ejection. 
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Figure 6. The SSE timestep at the moment when orbit- 
ing bodies are most likely to escape. Although short on the 
timescalc of Solar System evolution, this timestep is long rel- 
ative to a Solar rotation period, where additional variability 
which is not modeled here might occur. The mass loss over 
this timestep is treated as linear. 



3.6.2 Example Bounded Orbits for t]q < 0.5 

We consider here short orbits (ao = 3 x 10 3 AU) which are 
significantly disrupted by mass loss that takes place en- 
tirely within an orbital period (tjq = 0.3, see Fig. [1} and 
remain ellipses. Figure [7] illustrates the orbital evolution 
of 16 objects, all with different fo values, which all begin 
the N-body simulation with small-eccentricity (eo = 0.2) 
orbits. These objects may be treated as low-eccentricity, 
inner Oort Cloud bodies. Despite the small eo, the semi- 
major axes all increase by a factor of « 1.5 — 2.2, and the 
final values of e range from w 0.1 — 0.6. If this mass was 
lost adiabatically, then the semimajor axis would increase 
by a factor of m 1.5, the lower bound of the actual in- 
crease in a. Further, the fo = 240° object (brown curves) 
becomes nearly circularized at t « 12461.67 Myr; if mass 



loss ceased at this point, then the body would maintain 
this nearly circular orbit. 

Mass loss has a dramatic effect on the evolution 
of /. A system no longer is adiabatic when df/dt — 
()Veras et alJ l201ll ). At this point, / stops circulating 
and begins to librate. The lower left plot illustrates how 
the mass loss slows df/dt for each orbiting object. Al- 
though none of the orbiting objects achieve df/dt — 0, 
the fo = 90° object (black curves) comes closest. This is 
also the object whose a and e values are increased by the 
greatest amount. Alternatively, the highest ending value 
of df/dt is associated with the fo = 240° object, whose 
eccentricity decreases the most. In the linear mass loss ap- 
proximation, for $ 1, bodies closest to pericenter are 
most susceptible to escape. This plot helps demonstrate 
how no such correlation exists if the evolving regime is 
not runaway. 

Although the bottom-right plot might suggest that 
some planets appear to be leaving this system, all bodies 
will remain bound, orbiting the Solar white dwarf with 
their values of a and e attained at t = 12461.9 Myr, 
assuming no additional perturbations. 



3.6.3 Example Escape Orbits for r)@ ^ 0.5 

Here we consider longer orbits (ao = 1 X 10 4 AU) which 
are significantly disrupted by mass loss that takes place 
over several orbital periods (770 = 0.5, 0.6, see Fig. [1} and 
which can cause an orbiting object to escape. We consider 
highly eccentric orbits (eo = 0.8) and compare two types 
of escape in Fig. [8] the left panels showcase the fo = 0° 
(red) object achieving an orbit with e = 0.99872 and a > 
10 6 AU (two orders of magnitude higher than ao), while 
the right panels showcase the fo = 160° (pink) object 
achieving a hyperbolic orbit. The only initial parameter 
changed in the two plot columns was the value of t/q (0.5 
on the left, and 0.6 on the right). 

As evidenced by the plots in the upper two rows, the 
first strong burst of mass loss at t m 12.22 Gyr provides a 
stronger orbital kick for t]q = 0.6 Solar evolution. How- 
ever, as indicated by Fig. [2] the second, stronger burst of 
mass loss at t » 12.36 Gyr yields a lower value of a cr i t for 
t)q = 0.5. Regardless, this second strong burst of mass 
loss for both 7/0 values are roughly comparable, as the 
final range of a and e values are similar. 

The difference is enough, however, to cause escape in 
two different ways, one with an object which was at peri- 
center at the beginning of the simulation, and one with 
an object which was near-apocenter. This figure provides 
further confirmation of the difficulty with correlating fo 
with / for objects evolving in the transitional regime be- 
tween the adiabatic and runaway regimes. Both / vs t 
plots demonstrate the expected change of behavior dur- 
ing the second strong mass loss event. However, they fail 
to distinguish between both types of escape, and nei- 
ther technically reach df/dt = 0. However, for the left 
and right panels, respectively, df/dt = 0.068° /Myr and 
df/dt = 0.078° /Myr at t = 12461.9 Myr. 

The bottom panels show the spatial orbits of the ob- 
jects; the left panel illustrates all 16, and the right panel 
features only 4 (for added clarity), including the escap- 
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ing object. Due to the number of simulation outputs, 
lines are not connected between the outputted points. 
These plots explicitly illustrate how two bursts of mass 
loss from the Sun cause two distinct orbit expansions. 
Additionally, the apparent "noise" about the equilibrium 
orbits attained before the second mass loss event demon- 
strate the gradual expansion of the orbit due to minor 
mass loss events. Comparing the colours and positions of 
both plots shows almost no correlation with the identical 
initial conditions. The escaping pink body on the right 
panel never reaches a distance to the Sun which is closer 
than its initial pericenter value (2000 AU); the two dots 
for which x > 1.5 x 10 4 AU indicate the path along which 
the body escapes. 



3.6.4 Example Orbits for Rampant Escape 

We now consider objects in considerable danger of es- 
cape. These can reside at typical Oort Cloud distances, 
ao ~ 10 5 AU. The strong mass loss from the t)q =0.2 
evolutionary track (with maximum a ~ 7.5xlO- 6 M Q /yr; 
see Fig. [5} yields $ > 30 S> 1 and places these objects in 
the runaway regime. Here we determine if, unlike in the 
adiabatic regime, we can link the fate of an object with 
/o amidst nonlinear mass loss. 

Figure [9] confirms that we can. The results agree well 
with the linear theory. At the highest eccentricities sam- 
pled (eo = 0.95), the objects initially at apocenter (or- 
ange; /o = 180°) remain bound to the dying Sun, and, 
further, suffer the greatest eccentricity decrease. Addi- 
tionally, 4 other objects nearest to the apocenter of their 
orbits also remain bound. As the initial eccentricity of 
the objects is decreased (from 0.95 in the top panel to 
0.6 in the bottom panel), more bodies which are increas- 
ingly further away from their initial apocenter are allowed 
to survive. Also, as the initial eccentricity is decreased, 
more bodies which are progressively further away from 
their initial apocenter suffer a net eccentricity decrease. 
Note finally that the results are not symmetric about 
fo = 180°: objects approaching apocenter are slightly 
more protected than objects leaving apocenter. 



4 APPLICATION TO THE SCATTERED 
DISC, OORT CLOUD AND OTHER 
TRANS-NEPTUNIAN OBJECTS 

4.1 Overview 

The Solar System architecture beyond the Kuiper Belt is 
often divided into two distinct regions: the scattered disc 
and the Oort Cloud. The former is possibly an ancient, 
eroding remnant of Solar System formation, whereas the 
latter is thought to dynamically interact with the local 
Galactic environment and is continuously depleted and 
replenished. Both regions may supply each other with 
mass throughout the lifetime of the Solar System. 

However, the definitions of these regions are not stan- 
dardized. The boundary of the scattered disc is empirical 
and refers to objects with perihelia beyond Neptune (q > 
30 AU) and semimajor ax es beyond the 2:1 resonance 
with Neptune (a > 50 AU) (|Luu et al.ll 19971 : iGomes et al.l 



120081 ). Other papers have set a variety of numerical 
bounds based on their empirical readings of the semima- 
jor axis - eccentricity p hase space for object s beyo nd Nep- 
tune. In an early study, [Duncan fc Levisonl l| 19971 ) distin- 
guished the Kuiper Belt from a scattered disc of objects 
with q « 32 - 48 AU, e sC 0.8 and i < 50°. Alternatively, 
sometimes the scattered di sc is considered to be a sub- 
class of the Kuiper Belt fe.g. lGladman fc Chanll2006r ) and 
contains objects which satisfy the sim p ler constraint q > 
33 A U (|Tiscareno fc Malhotral 120031: lyolk fc Malhotral 
| 2008l) and e > 0.34 (|Volk fc Malhotral |200SD, whereas 



iKaib fc Quiiml (|2008r i give q ^ 40 AU and ISheppardl 
( 20ld ) claim q ~ 25 — 35 AU as th e criterion for in- 
clusion. Further, lLevison"et ID l|2004D demonstrate how 
the scattered disc population has changed with time. For 
example, at an age of 10 6 yrs, m ost scattered d i sc obj ects 
had q < 10 AU. More generally. ILevison et all j2008l ') de- 
fine the currently-observed scattered disc to be composed 
of objects with perihelia close enough to Neptune's orbit 
to become unstable over the (presumably main-sequence) 
age of the Solar System. However, some scattered disc 
objects may not be primordial, and instead represent a 
small transient population supplied by the Kuiper Belt. 
Also, because scattered disc objects interact with Nep- 
tune, some are depleted and diffuse into the Oort cloud. 

If the scattered disc is insufficiently replenished from 
the Oort Cloud and Kuiper Belt, then the scattered disc 
may be depleted relative to its current level by the end 
of the Sun's main sequence. Whatever disc objects which 
might remain will continue to interact with Neptune as 
their orbits move outward. If the orbits expand adiabat- 
ically, the scattered disc object will continue to interact 
with Neptune as during the main sequence. If the orbits 
expand non-adiabatically, then the scattered disc objects' 
eccentricity will change, and may no longer interact with 
Neptune. However, the semimajor axis of a scattered disc 
object which is high enough to cause non-adiabatic evo- 
lution is generally too high (a > 10 3 AU) to remain in 
that population. Therefore, the scattered disc population 
will be subject to adiabatic mass loss evolution, and con- 
tinue to_jnteract_by diffusing chaotically in semimajor 
axis (|Yabushitalll980l ). likely helping to repopulate the 
Oort Cloud. 

Beyond the scattered disc resides the Oort Cloud, 
which was originally postulated to contain ~ 10 11 comets 
of observable size with 5 x 10 4 AU < a < 1.5 x 10 5 AU, 
< e < 1 and isotropically-distributed inclinations (|Oortl 
1950). Some of these comets occasionally enter the in- 
ner Solar System and can achieve perihelia of a few AU. 
These intruders are thought to have been jostled inward 
by forces external to the Solar System. Many subsequent 
investigations since this seminal work have refined this 
basic picture for the outermost region of the Solar Sys- 
tem. 

The external perturbations ma y arise from passing 
stars and/or the local Galactic tide. iHeisler fc Tremaing 
(1986) demonstrated that tides likely dominate in- 
jection events into the inner Solar System, while 
passing stars continuously randomiz e the Oort Clou d 
comet orbital paramete r distribution l|Dvbczvriskill2002l ) . 
M orris fc Mulled ( 1986) demonstrate how Galactic forces 
vary the orbital angular momentum of an Oort Cloud 
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Figure 7. The orbital evolution of objects at o,q = 3 X 10 3 AU with eo = 0.2 from 12461.5 Myr. All of these bodies remain bound 
to the dying Sun, despite sometimes severe disruption to their orbits; the double cross in Fig. [2] corresponding to this initial pair 
(ao,??©) indicates that escape occurred for eo 0.8. The lower-right plot demonstrates that the disruption occurs within a single 
orbital period, and can stretch or contract orbits (upper-right plot) while widening them (upper-left plot). After this moment in 
time, the orbits assume osculating unchanging ellipses, and remain on these orbits as the Sun becomes a white dwarf. The lower-left 
plot indicates that all objects barely remain in the adiabatic regime because df(t)/dt ^ for all cases at all times, ensuring that / 
continues to circulate, albeit slowly. The orbit which is most severely disrupted (/o = 90°, black curve) comes closest to achieving 
df{t)/dt = 0. 



comet during a single orbit, which, hence, may be neither 
an ellipse nor closed. Therefore, if a typical comet's bound 
orbit significantly varies from Keplerian motion due 
to ext ernal perturbations, the analytics in IVeras et al.l 
(|201 lh cannot be applied to these objects without incor- 
porating models of those perturbations. 

However, if we do treat bound Oort Cloud cometary 
orbits as approximately elliptical, then we can estimate 
the effect of stellar mass loss on the cloud. One major 
refinement of Oort's original model is the division of the 
Oort Cloud into a n "inner" a nd "outer" population at 
a = 2 x 10 4 AU (|Hillsl Il98lh . This bifurcation point, 
which determines if a comet is observable only when it is 
part of a shower, is still widely utilized in modern sim- 
ulations and computations. Figure [2] indicates that this 
bifurcation value is a factor of 2 — 20 times as high as 
the Solar System's post main-sequence escape boundary. 
Hence, comets from both populations will be susceptible 
to escape due to post m ain-sequence mass lo ss. More- 
over, in a seminal study, iDuncan et all (|l987f ) approxi- 



mated t he inner e dge o f the Oort Cloud to be « 3000 AU, 
whereas iBrasserl (2008) use this same value to divide the 
Oort Cloud further into a third "innermost" region. High- 
eccentricity bodies with semimajor axes approximately 
equal to this potential second bifurcation point in Oort 
Cloud structure may be susceptible to escape during the 
Sun's post-main sequence evolution if r/Q < 0.6 (see Fig. 

Many modern Oort Cloud investigations focus on 
formation and/or the orbital properties of observed 
comets instead of topics more relevant to this study, such 
as the eccentricity distribution of the vast majority of 
the remaining (currently unobservable) comets and the 
future evolution of the Oort Cloud over the next severa l 
Gyr. However, som e authors (|Emel'Yanenko et al.ll2007l ; 
iBrasser et alj2oTol ) present plots which illustrate that the 
eccentricity distribution of Oort Cloud comets do span 
from to 1, with the minimum eccentricity typically in- 
creasing as a decreases. Therefore, the orbits studied in 
Sections 3.6.2 and 3.6.3 represent plausible families of 
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Figure 8. The orbital evolution of objects at ao = 1 X 10 4 AU with eo = 0.8 for r/Q = 0.5 (left panels) and r/Q = 0.6 (right panels) 
from 12.2 Gyr. In the left panels, the /o = 0° (red) body escapes when, after achieving e = 0.99872, its semimajor axis exceeds 
10 6 AU. Alternatively, in the right panels, the /o = 160° (pink) body escapes by achieving a hyperbolic orbit (e > 1, a < 0). 
The / vs t panels demonstrate that the true anomaly cannot distinguish between these two types of escape in this case, and that 
bodies may escape just before achieving df(t)/dt = 0. The pericentric vs. near-apocentric /o values of the escaping bodies, as well 
as well as a comparison of the orbit locations in the bottom two plots, indicate that the phase information is easily lost on the 
approach to df(t)/dt = 0. The bottom plots also demonstrate how the orbital architecture changes over several orbital periods, in 
contrast to Fig. [7] Note that the ejected planet in the bottom-right plot never comes within its initial pericenter value (2000 AU) 
of the Sun. 
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Figure 9. Eccentricity evolution in the runaway regime for 
typical Oort Cloud distances of 10 5 AU and t]q = 0.2. Objects 
initially closest to the apocenter of their orbit are the most 
likely to survive, and bodies at their apocenters experience the 
largest net drop in eccentricity. Objects further away from the 
apocenter may survive if their initial eccentricity is decreased 
(from 0.95 in the top panel to 0.6 in the bottom panel). 

"inner" Oort Cloud orbits, and the orbits presented in 
Section 3.6.4 represent a plausible family of "outer" Oort 
Cloud orbits. 

Multiple subpopulations of t he Oort Cloud help cat- 
egorize recent discoveries. Sedna (|Brown et al.|[2004h . the 
most distant Solar System object yet observed (at 90.32 
AU), highlights the failure of the traditional scattered 



disc and Oort Cloud to partition the Solar System be- 
yond the Kuiper Belt. With q « 76 AU and a » 480 AU, 
Sedna does not fit inside either population, and demon- 
strates that the orbital parameter phase space of trans- 
Neptunian objects is larger than previously thought. If 
Sedna maintains its orbit until the Sun turns off of the 
main sequence, then the object is guaranteed to remain 
bound (Fig. during post-main sequence mass loss, re- 
gardless of its position along its orbit. However, other 
objects a few hundred AU more distant might not sur- 
vive. Two other objects which de fy easy classification 
are 2006 SQ372 (|Kaib et al l l2009h . with q sa 24 AU 
and a « 796 AU , and 2000 00 67 (|Veillet et al.1 l200ll ; 
iMillis et al.ll2002r i. with q w 21 AU and a « 552 AU. The 
perihelions and aphelions of both objects imply that they 
can be classified as either scattered disc objects or Oort 
Cloud comets. Regardless, the short stability timescale 
(~ 200 Myr) predicted for 2006 SQ372 implies that it will 
not survive the Sun's main sequence phase. 2000 OOo7 
will remain bound (Fig. [2j to the post- main sequence Sun 
if the object survives the Sun's main sequence. 



4.2 Depletion Characteristics 

Although determining the fraction of the Oort Cloud 
which is depleted likely requires detailed modeling and 
extensive N-body simulations, we can provide some es- 
timates here. We perform an additional 3 sets of simu- 
lations, for 770 = 0.2, 0.5 and 0.8, and focus on the re- 
gion beyond 10 4 AU, where most Oort Cloud objects 
are thought to reside. These values of r)@ represent, in 
a sense, one nominal and two extreme Solar evolutionary 
tracks. For each we integrate with exactly the same 
conditions as described in Section 3.3, but now sample 
72 uniformly-spaced values of /o (for 0.5° resolution), 9 
uniformly-spaced values of eo in the range [0.1 — 0.9], and 
the following 6 values of a /AU= 1.0 x 10 4 , 2.5 x 10 4 , 5.0 x 
10 4 , 7.5 x 10 4 , 1.0 x 10 5 , 1.25 x 10 5 . 

The results of these simulations are presented in 
Figs. UOllllI Fig. [10] is a cartoon which illustrates the re- 
lation between instability and initial conditions, and Fig. 
[TJJplots the fraction of the Oort Cloud lost depending on 
the Solar evolutionary model adopted. The cartoon in- 
dicates unstable systems with red crosses, as a function 
of ao, eo and /o- The connection of these variables with 
instability helps describe the character of the mass loss. 
For example, the initial conditions which lead to insta- 
bility for the r/0 = 0.2 simulations (upper panel of Fig. 
1101) are patterned and focused towards the pericenter for 
higher ao and eo values, mimicking runaway regime be- 
haviour and implying large mass loss rates (see Fig. [SJ . 
In contrast, the initial conditions which lead to instabil- 
ity for the 770 = 0.5 simulations (middle panel of Fig. 1101) 
exhibit little order, implying that the mass loss rates are 
not large enough to reside fully in the runaway regime, 
but are large enough to allow for escape. The r}@ = 0.8 
simulations (lower panel of Fig. I10|) provide a third alter- 
native, showing a greater sensitivity to ao than to either 
eo or /o. The bottom two panels indicate the difficulty 
in predicting the prospects for escape amidst nonlinear 
mass loss without appealing to N-body simulations. 
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There are other features of Fig. [TUlworth noting. Not 
shown are objects with eo ^ 0.2, because all remain sta- 
ble. This implies that in no case does the Sun lose enough 
mass for a long enough period of time to eject distant 
circular bodies. In the upper panel, at least one object 
for every other (ao,eo) pair becomes unstable. Also, in 
this panel, asymmetric signatures about the pericenter 
indicate that of two bodies at equal angular separations 
from pericenter, the body approaching pericenter is more 
likely to escape. This asymmetry does n ot arise in the 
analytic linear theory (|Veras et al.l l201lT ). but may be 
guessed based on physical grounds, as a planet initially 
approaching pericenter is more likely to be ejected than 
a planet initially approaching apocenter. This tendency 
is also suggested in Fig. [9] with the top two curves in 
each panel. In the bottom-rightmost ellipse in the upper 
panel of Fig. 1101 we superimposed two brown diamonds 
at the locations of / cr it and 360° — / cr jt in order to deter- 
mine how closely the escape boundaries mimic the ana- 
lytic boundaries which ensure that a body's eccentricity 
must initially decrease. 



The middle and bottom panels of llOl indicate that the 
fraction of bodies which escape does not necessarily scale 
with ao; the bottom panel displays a dramatic double- 
peaked red cross-dominated structure at ao = 5 x 10 4 
AU and ao = 1 x 10 5 AU. One can glean understanding 
of those panels by noting three points: First, the orbital 
periods of objects in both panels is roughly equal to the 
mass loss timescales achieved for those Solar evolution- 
ary tracks, placing the objects in a transitional regime. 
Second, the rate of change o f eccentricity, whi ch is pro- 
portional to (e + cos/) (see IVeras et all 12011). can as- 
sume either positive or negative values and increase or 
decrease depending on the value of ao in this transitional 
regime. This result is not exclusive to nonlinear mass loss 
prescriptions; one can show that even for a linear mass 
loss rate, the inflection point of the eccentricity with re- 
spect to time contains multiple extrema in the transi- 
tional regime. Third, the moment when the Sun becomes 
a white dwarf and stops losing mass freezes an object's 
orbit, even if the orbit was approaching e —¥ 1. There- 
fore, if the rjQ = 0.8 post main-sequence mass loss lasted 
longer, the distribution of unstable systems in the bottom 
panel would no longer be bimodal. 



Fig. [TT]simply counts the unstable systems from Fig. 
1101 and both figures can be compared directly side-by- 
side. In the top panel of Fig. 1111 over 80% of the high- 
est eccentricity objects escape, and over half of all ob- 
jects with eo = 0.5 escape. The escape fractions become 
nearly flat for a^5x 10 4 AU. In contrast, in the middle 
panel, all eo curves except one maintain escape fractions 
which are < 20% for all values of ao that were sampled. 
In the bottom panel, the escape fractions spike at both 
a = 5 x 10 4 AU and a = 1 x 10 5 AU for the five high- 
est eccentricity curves. The results indicate how sensitive 
Oort Cloud depletion is to both the Oort Cloud model 
adopted and the Solar evolution model adopted. 



5 DISCUSSION 

5.1 Planetary Perturbations 

Perturbations from surviving planets which reside within 
a few tens of AU of the Sun will have no effect on the 
orbits described here, unless an external force is evoked, 
or one or more of the planets is ejected from the sys- 
tem and achieves a c l ose en counter on its way out. As 
shown in IVeras et all l|201ll l the pericenter of an orbit- 
ing isolated body experiencing isotropic mass loss cannot 
decrease. This holds true in both the adiabatic and run- 
away regimes, for any mass loss prescription. Therefore, 
mass loss cannot lower the pericenter of an Oort Cloud 
comet which was originally beyond Neptune's orbit. 

How the known surviving planets evolve them- 
selves under their mutual gravitational attraction amidst 
post-main sequence mass loss is less clear. With no 
mass loss, the giant planets are st able until th e 
Sun turns off of the main sequ ence l|Laskarl I2008T ); 
iKholshevnikov fc Kuznetsovl (|201ll ) indicate that Jupiter 
and Saturn need to be about 20 times more massive than 
their curr ent values in order t o suff er close encounters. 
However, iDuncan fc Lissauerl |l998T) demonstrate that 
with post main-sequence mass loss, giant planet evolu- 
tion may not be quiescent. In one linear mass-loss ap- 
proximation, they find that Uranus and Neptune's orbits 
will eventually cross. More generally, if the planets remain 
in orbits under ~ 100 AU and beyond the Sun's expand- 
ing envelope, they will evolve in the adiabatic regime, 
expanding their near-circular orbits all at the same rate. 
They will remain bound according to Figure [21 However, 
as multiple planets expand their orbits and maintain their 
relative separations, the critical sep arations at which the 
plane ts will remain stable will vary (|Debes fc Si gurdsson 
I2002T ). Also, as alluded to in Section 1, Jupiter and Sat- 
urn are predicted to largely maintain their current orbits 
until the Sun turns off of the main sequence. If this is 
true, then they will still be close to the 5:2 mean motion 
commensurability when significant stellar mass loss com- 
mences. How near-resonance behavior is linked to a po- 
tentially varying critical separation is not yet clear, and 
suggests that much is yet to be discovered about post 
main-sequence planetary systems with multiple planets. 

5.2 Beyond the Solar White Dwarf 

After the Sun has become a white dwarf, surviving ob- 
ject orbits will retain the orbital parameters achieved at 
the moment the Sun stopped losing mass. These objects 
will forever remain on these orbits unless subjected to 
additional perturbations. For the surviving planets, the 
primary perturbations will come from one another, which 
could lead to instab ilities on Gyr, or long er, timescales 
(see, e.g., Fig. 29 of IChatteriee et alH2008T ). For objects 
further afield, the dominant perturbations will be exter- 
nal, and come from, for example, Galactic tides (e.g., 
iTremainel Il993h Q and passing stars (e.g. IVeras et all 

1 which scales as M*(t)- 2 / 3 l lTremainj|l993l) and might har- 
bor an entirely different functional form from the present- 
day prescription due to the collision of the Milky Way and 
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Figure 10. Initial condition stability snapshot for r)Q = 0.2 
(upper panel). TjQ = 0.5 (middle panel) and r?Q = 0.8 (lower 
panel). Each scaled ellipse is for a given triplet (»?0, ao, eo) and 
graphically represent the initial fo of unstable systems (large 
red crosses) . The orange dot is a representation of the Sun, and 
the brown diamonds (at r)Q = 0.2 and the highest values of ao 
and eo) represent the values / cr i t = 180° — (1/2) cos -1 (7/9) ~ 
160° and 360° — / cr it ~ 200° , which, for runaway mass loss, 
bound the region in which the eccentricity must initially de- 
crease. For r)Q = 0.2, strong "runaway" mass loss causes ob- 
jects which are closest to pericenter to be lost more easily; the 
effect is enhanced for higher ao and eo - For ?7q = 0.5, the tran- 
sitional regime between adiabatic and runaway dominates, and 
objects appear to go unstable at a variety of nonconsecutivc 
values of fo. For tjq = 0.8, the instability appears to be highly 
sensitive to ao- 
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Figure 11. Fraction of objects which escape from the Sun 
for t]q = 0.2 (top panel), 0.5 (middle panel), and 0.8 (bot- 
tom panel) for the simulations described in Section 4.2 and 
shown in Fig. \W\ The e = 0.9,0.8,0.7,0.6,0.5,0.4, and 0.3 
curves are given by the colored curves with open triangles, 
open diamonds, open squares, open circles, downward filled 
triangles, upward filled triangles, and filled diamonds, respec- 
tively. Other curves do not deviate from 0%. Note the stark 
difference in all three distributions, demonstrating the sensi- 
tivity of Oort Cloud depletion to both the Solar evolutionary 
model and Oort Cloud model adopted. 
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2009). These perturbative sources will likely play a larger 
role post main-sequence than during the main sequence 
because of orbit expansion due to mass loss. These effects 
may eventually unbind the object. 

Define a cxt to be the minimum semimajor axis at 
which an external perturbation will eventually (at some 
arbitrary future time) remove an orbiting object. The 
brown dashed curve in Fig. indicates that the Sun will 
lose « 48% of its current mass. Therefore, an object cur- 
rently with a semimajor axis of ao 0.52a cx t evolving 
adiabatically during post-main sequence Solar mass loss 
will be eventually be subject to escape from these pertur- 
bations. However, if the object is in the runaway or tran- 
sitional regime, then ao will increase by a greater amount. 
For example, the red body in the left panel of Fig. [8] in- 
creases its semimajor axis by two orders of magnitude. 
In the idealized case of linear mass loss for a purely run- 
away object on a circular orbit at apocenter, Eq. (44) of 
IVeras et al.l (|201ll ) gives ao = 0.08a ex t assuming that the 
Sun loses 48% of its mass entirely non-adiabatically in 
the runaway regime. 



5.3 Exoplanet Analogues 

Although stellar evolutionary tracks are highly depen- 
dent on the zero-age main sequence mass and metallicity 
of a star, generally the maximum value of a increases as 
the stellar mass increases. Therefore, for stars of approx- 
imately Solar metallicity which are more massive than 
the Sun, a cr it is expected to be lower than for the So- 
lar System, and orbiting exobodies would be more prone 
to escape. The methodology in this work can be applied 
to any extrasolar system for which a stellar evolutionary 
track can be estimated. Given recent detections of ob- 
j ects which may b e orbiting parent sta r s at a > 10 3 AU 
llBeiar et all 120081 ; iLeggett et ail 120081 ; lLafreniere et al.l 
120111 ). this type of analysis may become increasingly rel- 
evant. 



6 CONCLUSION 

The Solar System's critical semimajor axis within which 
bodies are guaranteed to remain bound to the dying Sun 
during isotropic post-main sequence mass loss and in the 
absence of additional perturbations is ~ 10 3 AU - 10 4 
AU. A more precise value can be obtained for a given So- 
lar evolution model; this range encompasses many real- 
istic evolutionary tracks. The most important Solar evo- 
lutionary phase for dynamical ejection may be the RGB 
(tjq ^ 0.7), the thermally pulsing AGB (779 ^ 0.5) or the 
early AGB (0.5 < r/ Q < 0.7). Objects with a > a cr it, 
such as Oort Cloud comets, and those with a > a cr it, 
such as some scattered disc objects, may escape depend- 
ing on their orbital architectures and position along their 
orbits. Quantifying the fraction of the population of these 
objects which escape would require detailed modeling of 



Andromeda before the Sun turns off of the main sequence 
lICox fc Loebll2008t) . 



their (unknown) orbital parameters at the start of the 
Sun's post-main sequence lifetime. 
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